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An improved stale equation for the vicinity of the critical point is proposed. An analysis of the ex- 
perimental data on helium and xenon has been carried out in order to investigate the influence of the 
number of constants in the equation and the PpT ran^c on the critical constants T, and p, and on the 
critical exponents a, (3, y. and o\ No such influence has been detected. The model for the critical point, 
recently proposed by Widom, has been checked regarding its consequences for the rectilinear diameter. 
No definite confirmation but indications for its correctness have been found. 
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1. Introduction ru ,. c dp m , . ,. ... 

1 he divergence <>t ~~tf , where p m is the rectilinear 
<II 
Recently an equation of state for the vicinity of the diameter, is a consequence of a model for the critical 
critical point has been proposed by Verbeke et al. [1 J. 1 state, which has recently been proposed by Widom 
This equation was applied to the data for xenon of and Rowlinson [5]. This divergence is given by another 
Habgood and Schneider [3]. In a subsequent paper the power law, 
equation was applied to data for methane of Jansoone 

et al. [2]. It is the purpose of this paper to report im- dp m ^ _ _ aj 

provements in the equation which have been made dT 

with a\ = a (as derived by Widom). The same model 

2. Formulation of the Modified State Equation yields a similar divergence of^, where p, is the line 

ci • i • • l .• of symmetry in the CAP, T) surface where C r is the 

several questions concerning the original equation •/» i_ i a ■• r 

remained to be answered. »l* l ' lfic heat at con8tanl volume - A divergence of 

The temperature derivative, ^pf, where P, is the ^ was incorporated in our original eq (1) in order to 

i rp . , -ci- j provide a similar divergence in all terms of that equa- 

vapor pressure and I is the temperature, is believed J. • \ • *u * *• i 

*V ^ r tion; since A in that equation was not necessarily 

& equal to 1, the corresponding ct\ was restricted to 

^~{T,--T)-\ (1) a,= l-j8. 

close to the critical point. The value of a 2 must, accord- It was found, however, that a value of k different from 

ing to the weak inequality of Griffiths which is based on 1 does introduce a discontinuity across the critical 

thermodynamics, be less than or equal to a-\- /3, where isochore. 

/3 is the exponent of the power law for the difference In order to test the conjectures on ol\ and a>> made by 

of the coexisting densities and a is the exponent of the M. S. Green et al. and Widom respectively, our 

power law for C v . It has been proved [7] for an equation original equation [lj was modified as follows: 

of state of the type proposed by M. S. Green et al. [4], 

that a must be equal to a 2 . The present equation can P= P*+ B(p) (T—T s ) + Ap(p — p t ) [(T 

also, but must not necessarily, imply the latter equality. + T c — 2T s ) i+fX — (T c — T 8 ) 1+M ] (2) 

with 
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B(p)=P ll + Bi(p-p m )\p-p m \^ + B,. . ., 

T, = T c -[C\p-p m \yie 



and 



pt = pc + pto(T-T s y-«>. 



(4) 
(5) 
(6) 

(7) 



It can be shown from eq (2) that the exponent, a, of 
the power law for C v can be expressed in terms of 
/jl and /3 as follows: 

a=l-(2j3+/x). 

Several questions can now be formulated. 

1. Is a = a\ = 0:2? 

2. Is p m \ equal to — p (i) ? 

3. Is there any PpT range-dependence for the 

exponents /3, y, 5 for the following power 
laws p e — p g ~ (Y C — T)P (the indices / and 
g refer to liquid and gas), 



P \dp/T 



\Tr-T\~ 



(at 



saturation or at the critical isochore), 

\p-Pc\t=t ( t \P-Pc\s- 
4. How many terms are relevant in equation (6)? 

Answers to these questions have been obtained by 
applying eq (2) to two substances for which excellent 
data is available, i.e., the data on helium by Roach [6] 
and the data on xenon by Habgood and Schneider 
[3], As far as the latter data are concerned, a rather 
large difference was detected previously between the 
measurements on isochores and on isotherms [1|. The 
error (for version 2) with the isotherm data included 
amounted to 0.00311 atmospheres versus 0.00227 
without. In the subsequent analysis only measurements 
on isochores have been considered. 

3. Analysis and Results 

Before estimating the PpT range dependence of 
the critical point parameters, the ranges have to be 
defined. The choice of the ranges will always be more 
or less arbitrary. In this case the ranges are defined 
by confining the variables density, temperature and 



pressure within a selected percentage around the 
cirtical point parameters. The method of defining 
the various ranges is illustrated in figure 1. The actual 
ranges indicated by the cross hatching in figure 1 
are given in table 1. A single significant parameter, /V, 
the number of data points in the range, results from 
this limiting procedure; N will be taken as a quanti- 
tative measure of the range. In table 1 the boundaries 
are defined for each range and for each substance, 
and the corresponding number of data points in the 
range is given. Range No. 1 for helium extends farther 
than any other range, and for this range the number of 
coefficients, B, in eq (6) was varied in order to investi- 
gate their relevance. 




FIGURE 1. Schematic P— T diagram showing boundaries of a 
"range''' around the critical point. 

The variation of the error and the apparent values 
of the scaling exponents a, /3, y, and 5 as a function 
of the number of coefficients in eq (6) is shown in 
figure 2. Z?o and B\ were found to provide a satisfac- 
tory description for data of the given accuracy within 
the given range. Beyond these two constants, no 
dependence on the number of Z?'s was found for a 
and y and only a slight dependence for /3 and 8. 
Consequently in further calculations all other B-s 
have been assumed to be zero. Equation (2) can be 
applied on the data in two alternative versions; version 
1 where ai, a-2, p»n and pto are independent adjustable 



Table 1. Boundaries of the different PpT ranges for xenon and helium and 
the number of data points left in each range 



No. of 
range 



/ lPmax-Pr| \ 



(IW51) 



(^) 



Number of data 
points, TV 



Xenon 



Helium 



1 
2 
3 
4 



0.50 
.35 
.30 
.10 



0.05 
.035 
.03 
.01 



0.08 


173 


.05 




.05 


142 


.02 


114 



414 
295 



208 



45 
44 
43 
4.2 

4 1 

40 

39 
i 

Fk 



fi 



fr~ 



.$— 



^ - -~* 



HELIUM 
(Range I .Version 2) 



*— — — _ 



Number of constants B 



URE 2. a, /3. y (//?»/ & and their dependence on the number of 
Wsin eq{2). 



parameters and version 2, where a\ = a-i= 1 — (2/3 + jjl) 
and pm\= — pm. The two versions are compared tor 
helium and for xenon in table 2. The difference in 
the quality of fit between the two versions is not sub- 
stantial, and therefore further speculations on the 
ran^e dependence of the critical parameters have 
been carried out with version 2. The dependence of 
the critical exponents a, o. y and (3 on the ran<2;e is 
illustrated in figure 3. The constants for range 1 and 
for version 2 of equation (2) are given lor helium in 
table 3, and for xenon in table 4. 

The estimation of the parameters has been carried 
out through an iterative procedure. Initial constants 
have been chosen and the equation has been linear- 
ized with respect to the deviations of the constants. 
These deviations were then calculated by the well 
known least-squares methods for the linear case. 

The error analysis is carried out on the basis of 
the correlation matrix for the linear case. There may 
be some doubt whether this is entirely justified in 
the case of a nonlinear equation. 



Table 2. Comparison of the relevant parameters for the two versions of eq (2) and for helium and xenon 







Helium 


Xenon 


Ranjje No. 


Parameter 






















Version 1 


Version 2 


Version 1 


Version 2 


1 


8 


4.28 ±0.09 


4.16 ±0.06 


4.45 ±0.14 


4.39 ±0.101 




y 


1.128 ± .004 


1.119 ± .003 


1.232 ± .012 


1.191 ± .006 




p 


.343 ± .003 


.353 ± .002 


.356 ± .005 


.351 ± .005 




a 


.184 ± .008 


.173 ± .006 


.054 ± .020 


.108 ± .014 




«i 


.270 ± .045 


.173 ± .006 


.287 ± .116 


.108 ± .014 




ot-i 


.400 ± .026 


.173 ± .006 


.470 ± .045 


.108 ± .014 




Pt» 


-.00100 ± .00015 


-. 00132 ± .00015 


-. 0000417 ± .0000131 


-.00001 45 ± .000004 




p,m 


+ .00208 ± .00022 


-p» 


+ .0000453+ .0000186 


-pr«» 




* 

/ SAP, 

V/V-/VC atm 


0.00072 


0.00074 


0.00224 


0.00227 


2(He) and 3(Xe) 


o 


4.24 ±0.24 


4.13 ±0.08 


3.73 ±0.76 


4.72 ±0.39 




y 


1.137 ± .012 


1.138 ± .009 ' 


1.178 ± .016 


1.208 ± .012 




(i 


.351 ± .009 


.363 ± .006 


.430 ± .052 


.324 ± .012 




a 


.161 ± .023 


.135 ± .015 


-.039 ± .106 


.142 ± .033 




a, 


.277 ± .286 


.135 ± .015 


-.149 ± .389 


.142 ± .033 




OL> 


.338 ± .124 


.135 ± .015 


.060 ±1.591 


.142 ± .033 




p ( u 


-.0003 14 ± .000439 


+ .00013+ .0004 


-. 0000105 ± .0000070 


-. 000025 1± .0000059 




Pm x 


+ .00168 ± .00107 


-p» 


+ .000228 ± .000057 


-p.. 




* 

/ SAP,- 

V/v-zvc 


0.00059 


0.00059 


0.00198 


0.00203 


4 


5 

y 

a 
ot\ 

« 2 

Pn> 

pm\ 








5.73 ±8.19 
1.284 ± .074 
.271 ± .129 
.172 ± .283 
.172 ± .283 
.172 ± .283 
+ .0000171+ .0000425 
— p t » 




* 












/ SAP,- 

V/v-/vc 


0.00209 



*N is the number of data points; NC is the number of eonstants. 



209 





Table 3 


Constants for eq (2) applied to helium data [6] a 




Constants of the equation 


The standard deviation is 0.00074 




Constant 


Value 


Error 




A 


0.32043251E 04 


0.10131254E 03 




P„ 


.17091894E 01 


.19034584E-02 






B l 


.33758179E 03 


.12381564 02 






Ps2 


.72719432E 00 


.34948345E-01 


Maximum error in pressure is —0.0025 


atm. 


P.sO 


-,66253524E 01 


.98947437E 02 


At T= 5.26370 and 




pto 


-.13174329E-02 


.15183463E-03 


P = 2.40487 atm. 




2-ol 2 


= 1 -h 2)8 + /x 








c 


0.74227144E 02 


.45068578E 00 






Pc 


.17232286E-01 


.52370344E-05 






pm\ 


-Pto 








T c 


.51909878E 01 


.27286947E-03 






P 


.35366522E 00 


.23593233E-02 






P 


.11915742E 00 


.33542782E-02 






1-a, 


= 2/3 + m 








Derived Constants 




Constant 


Value 


Error 




a 


0.17351214E 00 


0.60508733E-02 




y 


.11191574E 01 


.33542782E-02 






8 


.41644543E 01 


.62999387E-01 







a The density range is 50, the temperature range is 5 and the pressure range is 8 percent around the critical point. The number of data 
points is 414. 



Table 4. Constants for eq (2) applied to xenon data \3] i] 



Constants of the equation 


The standard deviation is 0.00227. 




Constant 


Value 


Error 




A 


0.94959356E 04 


0.21903837E 03 




Psi 


.11860651E 01 


.83952672E-03 






£, 


.13935997E 04 


.53686691 E 02 






Ps2 


.15611421E-01 


.23977382E-02 


Maximum error in pressure is 0.0067 


atm. 


P,o 


-.28601355E 03 


.24339522E 00 


At r=289.13990Kand 




Pto 


-.14547437E-04 


.39990010E-05 


P = 56.3430 atm. 




2-a 2 


= l + 2^ + /x 








c 


0.47872048E 03 


.12869899E 01 






Pc 


.84971537E-02 


.39913666E-05 






pm\ 


— pto 








T c 


.28976481E 03 


.31386983E-02 









.350500350 00 


.46079306E-02 






P 


.19110181E 00 


.63731809E-02 






I -a, 


= 20+/* 








Derived Constants 




Constant 


Value 


Error 




a 


0.10789750E 00 


0.13966633E-01 




y 


.11911018E 01 


.63731809E-02 






8 


.43982899E 01 


.10541327E 00 







a The density range is 50, the temperature range is 5 and the pressure range is 8 percent around the critical point, the number of data 
points is 165. 
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tlGURE 3. a. /3. y and 8 for helium and xenon and their dependence 
on the number of data points (error bars are indicated by the cross- 
hatched strips). 

4. Conclusions 

Only qualitative conclusions can be drawn from 
table 2. The values obtained for a x and a 2 , indeed, are 
consistent with the divergences as claimed by Widom 
and by Green et al. The values obtained, however, are 
larger than a. Furthermore the equality Pmi = —pto 
seems to check out qualitatively. It should be em- 
phasized however that the accuracy of the data, 
although being within the limits of the present state 
of the art, does not allow definite conclusions. The 
small differences obtained for a, /3, y, and 8 in versions 
1 and 2 and by other authors [7] could be explained by 
the fact that equations of this type can be very sug- 
gestive with respect to these values. 

It should also be noted that the value for a obtained 
for helium is much higher than the values obtained in 
the past by other authors [7]. 

In a more recent analysis, however, higher values 
have been obtained [8];' other recent experiments on 
helium also suggest higher values for a (e.g.. a = 
0.14 [9]). 



It can be seen from the error bars in figure 3 that 
helium is the more likely candidate for interpretation 
of the resulting exponents: 

The agreement between the exponents for the differ- 
ent ranges and for the different versions of the state 
equation for helium is almost within the error bars. 

Noteworthy is the fact that a x tends to be higher 
(~0.3±0.1) than a and that a 2 seems to be even 
higher (~ 0.45) for version 1 and for both helium and 
xenon. 

We do not exclude the possibility that these values 
are suggested by version 1 but it is noteworthy that an 
a 2 equal to 0.45 is still lower thana + /3 (0.15 + 0.355). 
This value (a 2 = 0.45) however, would imply a thermo- 
dynamic potential divergent at the critical point, a 
possibility which is doubted by many. 

A weak point in eq (2) is the nonanalyticity intro- 
duced by the second term on the right side of eq (6) 

resulting in a divergent (— -) for p equal to p m . This 

\op /t 

nonanalyticity can be avoided but we doubt that its 
influence on the critical parameters is important. 
More difficult to avoid is the nonanalyticity intro- 
duced by the third term of eq (2). This nonanalyticity 
disappears when y//3 is an integer, say yl/3 = 3. for 
instance. This effect may have had the odd result of 
producing rather high values for a. 



We are indebted to J. M. H. Levelt Sengers for 
stimulating discussions on this subject, and to V. 
Jansoone who suggested the error analysis. Most of 
of the calculations were performed on the CDC3800 
computer, the use of which was granted by the Cryo- 
genics Division of the National Bureau of Standards. 
Boulder, Colorado 80302. We are also indebted to the 
"Rekencentrum" of the University of Leuven. Belgium 
where part of the computations were performed. 
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